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We report the observation and quantitative characterization of driven and spontaneous oscilla- 
tions of quantum entanglement, as measured by concurrence, in a bipartite system consisting of a 
macroscopic Josephson phase qubit coupled to a microscopic two-level system. The data clearly 
show the behavior of entanglement dynamics such as sudden death and revival, and the effect of 
decoherence and ac driving on entanglement. 
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I. INTRODUCTION 

Entanglement is a unique property manifesting quan- 
tum correlation of multiparticle quantum systems that 
has no classical counterpart. It has been one of the most 
fascinating and nonintuitive concepts of quantum me- 
chanics and has stimulated extensive debate [H Re- 
cently, interest in entanglement has intensified since it 
is considered as one of the key resources for quantum 
information processing [31 0] and as a consequence a va- 
riety of properties of entanglement have been discovered 
[SHZ]- Nevertheless, many fundamental questions about 
entanglement remain open, including the entanglement 
of autonomous open quantum systems, the effect of ex- 
ternal driving on entanglement, and the mechanism of 
damped entanglement oscillation (DEO), entanglement 
sudden death (ESD) and ESD revival (ESDR). Another 
important issue in the experimental study of entangle- 
ment dynamics is to find simple methods to measure en- 
tanglement. 

Entanglement can exist not only in microscopic but 
also in macroscopic systems such as Josephson phase 
qubits ( JPQs) [SI [H| , which are basically current or flux 
biased Josephson tunnel junctions having Josephson cou- 
pling energy Ej much greater than charging energy Ec- 
JPQs are essentially manufacturable atoms whose Hamil- 
tonians can be custom designed and realized with in- 
tegrated circuit fabrication technology [T0Ul2j . This 
unique property makes the JPQ a good test-bed for 
studying fundamental issues in quantum mechanics and a 
promising candidate for implementing quantum informa- 
tion processing. In our experiment reported here a flux 
biased JPQ, which is a radio frequency superconduct- 
ing quantum interference device consisting of a super- 



conducting loop of inductance L « 770 pH interrupted 
by a 4.8 /im^ Josephson junction of capacitance C w 240 
fF and critical current Ic ~ 1.4 /xA, is used as shown 
in Fig. 1(a). The two lowest levels in the upper well 
of the strongly tilted double well potential form the two 
computational basis states |0) and |1). The energy level 
spacing, cjiq (for convenience we set H = h/2Tr = 1 where 
h is Planck's constant) , between the two basis states can 
be continuously tuned by varying the amount of mag- 
netic flux inductively coupled to the superconducting 
loop. Although atomic size defects in tunnel barrier of 
the Josephson junction, which are essentially microscopic 
two-level systems (TLSs), are considered as one of the 
major sources of energy relaxation and decoherence in 
JQP [131 m] , they also have the potential to be utilized 
as a useful resource for quantum information process- 
ing [TSHIS]. In this work, we used a TLS coupled to a 
JPQ to investigate the dynamics of entanglement of the 
coupled bipartite system. We show that for a bipartite 
system described by Hamiltonian ([T]) the degree of entan- 
glement, quantified by "concurrence" (20j, in both driven 
and free evolution states can be obtained by measuring 
the state of JPQ alone, which is quite different from the 
conventional tomography method in studying the entan- 
glement. The results clearly show that resonant ac drive 
and always-on qubit-TLS interaction together generates 
entanglement oscillations and that in the subsequent free 
evolution the system may exhibit DEO, ESD, and ESDR 
depending on the specific bipartite states at the time 
of turning off the ac drive. A comparison between the 
experimental results obtained with finite relaxation and 
decoherence and the predictions based on analytical solu- 
tion of the corresponding pure states (i.e., relaxation and 
decoherence free) clearly shows that not only the system- 
environment interaction but also the bipartite system's 
initial states play important roles in determining the dy- 
namic behavior of this type of open quantum systems. 
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II. EXPERIMENT 

Fig. 1(b) shows the schematics of the circuitry used. 
To distinguish qubit states from those of TLS, we denote 
the ground and excited states of the qubit (TLS) as |0) 
(|g)) and |1) (|e)), respectively. The coupled qubit-TLS 
in a weak microwave field is described by a four-level cou- 
pled bipartite system whose effective Hamiltonian, in the 
basis of the four product states {\0g) , \ lg) , |0e) , |le)}, is 
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where the diagonal elements are the energies of corre- 
sponding product states, flm is the Rabi frequency of the 
JPQ, g characterizes the qubit-TLS coupling strength, 
and to is the frequency of the microwave field. Notice 
that at the degeneracy point one has Eig = i?oe-The 
system's parameters in Hamiltonian ([I]) were determined 
from spectroscopy {2g/2TT = 76.0 ±1.2 MHz), Rabi oscil- 
lation (Sl^/2iT = 63.0 MHz), pump-probe {Tf « 61 ns), 
and pump-SWAP-delay-SWAP-probe (Tf ^'^ « 146 ns) 
experiments (17j as shown in Fig. 2. 

To measure dynamics of the coupled system, we first 
prepare the system in the ground state \0g) at i < 0, 
followed by applying a resonant microwave pulse of width 
imw at t = to coherently transfer the system from |0g) 
to other states through qubit-microwave coupling and 
qubit-TLS coupling g. After the pulse is terminated at 
t ~ tmw, the probability of finding the qubit in the state 
|1), Pi — Pig -\- Pie, is measured after a time thcc is 
elapsed from the end of the microwave pulse as shown 
in Fig. 1(c). The procedure is repeated for different 
imw and ifroc to obtain Pi as a function of tmw and free 
evolution time tfroc- 



III. RESULTS AND DISCUSSION 

To facilitate a direct comparison to analytical result in 
our experiment the qubit-TLS was biased at the center 
of the anticrossing where Eig — Eoe- Fig. 3(a) shows the 
complete set of Pi(imw, ^frcc) data measured. To clearly 
illustrate the effect of on vacuum Rabi oscillation 
of Pi at t > tmw, data taken with iniw = 7.0, 13.5,20.5, 
and 27.0 ns are plotted as symbols in Fig. 3(b). Notice 
that while the initial phase of Pi oscillation depends on 
^mw the frequency is independent of t„^^ and its value, 
76 MHz, agrees very well with the size of the splitting 
2g obtained from the spectroscopic measurement. Fur- 
thermore, the data with t^yj = 13.5 ns and imw = 27.0 
ns show very little oscillation while in a stark contrast 
the data with tniw = 7.0 ns and t^iw = 20.5 ns oscillate 
with much larger amplitudes. Their oscillations are out 
of phase with each other indicating the importance of 
initial phase of vacuum Rabi oscillation in determining 
concurrence as discussed below. 



To understand the dynamics of this coupled bipartite 
quantum system, we first discuss the ideal situation of 
pure state evolution by solving the problem analytically. 
In rotating frame Hamiltonian (1) is transformed to: 
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where 0(t„„ — t) is the heaviside step function with 
9{tmw-t) = 1 for i < tmw and 0{t^^-t) = for t > t^^. 
It is obvious from Hamiltonian ([2| that the dynamics 
occurs only in the subspace {llg), |0e)} for t > ^mw 
Writing the wavefunction of the system in the form 
\^{t)) = CQg{t)\Qg) + cig{t)\lg) + coe(i)|Oe) + cie(i)|le), 
we obtain probability amplitudes by solving the time de- 
pendent Schrodinger equation directly. For t ^ tmw) we 
have 
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with fl^ — fig + g, — Qs ^ g, and ilg = + 9^- 

Notice that the quantity measured directly in our exper- 
iment, Pi{t) = |cig(t)p -I- |cie(t)p, undergoes anomalous 
Rabi oscillation which in general contains all three fre- 
quency components ri+, and fig- 

When the system is driven by a resonant microwave 
field P2{t < tmv/), being the probability of finding the 
system in the subspace spanned by \lg) and |0e), under- 
goes sinusoidal oscillation: 
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which shows that amplitude of the oscillation depends 
only on the ratio g/Q^. In our experiment, g/flm — 
0.60 yielding P2 max — 0.74 which agrees well with the 
experiment. 

For t > tjiiw, ifrce = t — ^mw, it is Straightforward to 
show 
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In this case, the probabiUty of finding the qubit in sate 
|1) can be expressed as 

A(i >tmw) ^ |cio|^ + |cig|^ 

= -Ple(imw) + 2-P2(imw)(l + C0s{2gtf,cc + 5^mw)), 

(5) 

with Pie(tmw) = |cie(imw)P being the population 

of |le) and P2(imw) = Plg{tmw) + Poe{tmw) = 

{n^+n^)'^ sijj2 Oiinw being the probabihty of finding the 
system in the subspace spanned by \lg) and |0e) at t = 
imw which remains constant for t > tmw Eq- ^ shows 
that after microwave is turned off the system undergoes 
vacuum Rabi osciUation caused by the interaction be- 
tween \lg) and |0e). The angular frequency, depth (peak- 
to-peak), initial phase, and bottom envelope of the oscil- 
lation are 2g, P2(^mw), gtmw, and Pie(imw), respectively. 
In addition, because Pog(imw) + -P2(*mw) + -Pie(imw) = 1 
the difference between unity and the top envelope of 
Pi{tmw,thcc) is just Pog{t 

Interestingly, when ilstmw = ^mr, P2(imw) = 0, i.e., 
vacuum Rabi oscillation vanishes and one has 



(6) 



Pl{tmw,t[rec) = sin' 



which is independent of ifroo- In addition, if f2_(-t,„„ ~ 
{2k + l)7r, one has Pi ~ 1, which corresponds to the 
population mainly occupying |le); if fi+imw — 2fc7r , one 
has Pi ~ 0, which corresponds to the population mainly 
occupying \0g). These two cases correspond to tniw — 
13.5 ns and 27.0 ns, respectively, as shown in Fig. 3(b) 
where the exponential decay is due to energy relaxation. 

Although the basic physics is the same, because of de- 
coherence and relaxation the experimental system inves- 
tigated here cannot be represented by pure states but 
mixed states described by the bipartite density operator 
p (t) . The diagonal matrix elements Pmm and off-diagonal 
matrix elements pmn {fn ^ n) represent the occupation 
probability of the state \m) and coherence between the 
states \m) and respectively. To simulate the system's 
dynamics, we solve the master equation [3T1 
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where, L^^i j^i'^i' \Hmm'^n'n Hn'n^mm'\ t Hrnm' 

is matrix element of the system's Hamiltonian, and 
Rmn,m'n' IS the damping rate matrix element whose value 
is proportional to the energy relaxation rate [2T] . Eq. Q 
is numerically integrated to obtain p (t) . The results is 
shown in Fig. 3(c) as well as the solid lines in Fig. 2 
(resonant ac drive) and Fig. 3(b) (free evolution). It 
can be seen clearly from Fig. 2 that unlike the normal 
sinusoidal Rabi oscillation observed in the region of large 
qubit-TLS detuning, at the degeneracy point the oscil- 
lation of Pi{t < tmw) is clearly non-sinusoidal due to 



more complicated dynamics of the driven four-level sys- 
tem [23ti26j . In Fig. 3(d) we also show the calculated 
Pi based on the analytical solution Eq. ([5| by treating 
the effect of energy relaxation phenomenologically. No- 
tice that agreement between the experimental, numeri- 
cal, and analytical results in the entire range of driven 
and free evolution is very good confirming quantitative 
understanding of the system's dynamics. 

To study the dynamics of qubit-TLS entanglement, we 
examine how concurrence, denoted as Cp |20j for mixed 
states, evolves with time. In Fig. 4(a) Cp derived from 
the measured Pi is shown. Of particular interest is that 
Cp is observed to undergo damped oscillation in both of 
the driven and free run (i.e., autonomous) parts of the 
evolution. Fig. 4(c) shows the Rabi-like oscillation when 
the system is driven by the resonant microwave field. The 
corresponding Cp (the dashed line) oscillates and under- 
goes sudden death and revival repeatedly. The extrema 
of Cp are correlated strongly with the distinctive "shoul- 
der" feature in Pi. Fig. 4(b) shows Cp in the free run 
part of the system's evolution with t^w = 7.0, 13.5,20.5, 
and 27.0 ns, respectively. When i,,™ = 7.0 ns and 
taiw — 20.5 ns, the qubit and TLS are mostly in {\lg) 
and |0e)}. The coupling between these two basis states 
via the g{\lg){Oe\ + \Oe){lg\) term of the system Hamilto- 
nian thus leads to time-dependent entanglement causing 
Cp to oscillate. In this case, the entanglement dynamics 
clearly exhibits the phenomena of no ESD (NESD) and 
ESDR [27 , respectively. In contrast, when i^w = 27.0 
ns (^mw = 13.5 ns), the system is mostly in \0g) (|le)), 
which is decoupled from all other three basis states after 
microwave was turned off (neglecting energy relaxation) . 
This results in a weak entanglement characterized by a 
small overall value of concurrence Cp at the start of free 
evolution and, when decoherence is taken in account, the 
entanglement sudden death. 

Further insights into this type of system's entangle- 
ment dynamics and the effects of environment can be 
obtained by examining time-dependent concurrence of 
the corresponding pure state system and comparing it 
with the experimental result. For the corresponding pure 
states, concurrence C — 2|cogCie — cigCoel ^28J. From Eq. 
([3]), we obtain 
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Eq. ([s]) shows that when driven by a resonant microwave 
field entanglement oscillation is rather complex which 
has three frequency components f2+, fl-, and g. In the 
free evolution stage the so called preconcurrence |20j . 
which is the quantity inside the absolute value sign of 
Eq. ([8|, undergoes sinusoidal oscillation with amplitude 
P2itmw), which can be obtained directly from measured 
Pi(imw, ^frcc) according to Eq. ([s]), and a vertical offset 
|/(r2in, i,„„)|. Hence, concurrence exhibits the "high- 
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low" oscillation shown in Fig. 4(b) unless offset of pre- 
concurrence is zero. By varying and measuring the 
subsequent vacuum Rabi oscillation one can trace time 
evolution of concurrence of the driven as well as the au- 
tonomous stage of evolution. The result also shows that 
measuring the state of phase qubit alone is sufficient to 
gain all information about entanglement dynamics of this 
bipartite system. 

In particular, when Jlg^mw = 2n7r one has P2(imw) = 
as discussed above and thus only \0g) and |le) contribute 
to concurrence: 

C = 2|cogCie| = 2^Pi{l-Pi) = I singi„iw|, (9) 

which is independent of iftcc- In addition, when g^mw — 
(2fc + l)7r and gt^w — 2fc7r, where k is an integer, we have 
C <C 1 [notice the different vertical scales in Fig. 4(b)] 
because either Pi ~ 1 or Pi ~ 0. 

According to Eq. (|8]), for pure states the amplitude 
of concurrence oscillation remains constant for all t and 
C reaches zero at a discrete set of times only as shown 
in the insets of Fig. 4(b). In contrast, the experimental 
result displays a variety of interesting behaviors includ- 
ing DEO, ESD, and ESDR illustrated in Fig. 4(b) (from 
top to bottom). Because the main difference between the 
measured qubit-TLS system and the corresponding hy- 
pothetical pure state system is that the former is an open 
quantum system interacting with its environment (e.g., 
energy relaxation and decoherence) while the latter is 
isolated, our result shows that environment is the dom- 
inant mechanism of the observed complex entanglement 
dynamics [29j. 



IV. CONCLUSION 

In summary, we have experimentally demonstrated 
that the coupled qubit-TLS system is a test bed for quan- 
titatively studying the dynamics of bipartite entangle- 
ment. The measured time evolutions of this bipartite 
system, either when driven by a resonant microwave field 
or in free evolution, agree very well with the analytical 
and numerical solutions. Our results demonstrate that 
in situations similar to those described here one not only 
can quantify entanglement via concurrence by measuring 
the state of one constituent only, but also be able to con- 
trol the dynamics of the entanglement by adjusting the 
interaction time between the qubit and the ac resonant 
driving field. A comparison between the temporal evolu- 
tions of concurrence of the open and the corresponding 
isolated systems indicates that for the bipartite system 
studied here the entanglement oscillation and revival are 
originated from the qubit-TLS coupling while the entan- 
glement decay and sudden death are due to the coupling 
to the environment. 
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FIG. 2: (Color online) Spectroscopy of the coupled qubit and 
TLS versus flux bias. The splitting due to the coupling be- 
tween the qubit and TLS is 2g/2-R = 76.0 ±1.2 MHz at 
/ = 15.722 GHz. The inserts are Rabi oscillation, Ti, Rabi 
beating and vacuum Rabi oscillation, separately (from top to 
bottom). 



FIG. 1: (Color online) (a) SEM of the sample which is fab- 
ricated by Al/AlOx/Al trilayers. (b) Schematic of the qubit 
circuitry. Josephson junctions are denoted by the X symbols, 
(c) A time profile of manipulation and measurement. The 
corresponding potential energy landscape is also shown. 
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FIG. 3: (Color online) (a) Experimentally measured Pi ver- 
sus tniw and tfrcc- (b) Pi oscillation with typical tm„. The 
symbols are the experimental measured Pi and the lines are 
numerical result, (c) Numerical result of Pi versus t^w and 
tfroc- (d) Analytical result of Pi versus tmw and tfroo- 
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FIG. 4: (Color online) (a) Concurrence as a function of tmw 
and tfroc extracted from the experimental data, (b) Concur- 
rence with tniw=7.0, 13.5, 20.5, and 27.0 ns, respectively. An- 
alytical results of concurrence without decoherence are shown 
in the insets. Notice that the different vertical scales are used, 
(c) Pi and concurrence as a function of tmw Open circles are 
experimental data and the solid line is the numerical result of 
Pi. The dashed line is the corresponding concurrence. 



